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6.2 Micro-scale basis of seepage flow, 
Theory of homogenization 

The phenomenological law of Darcy describing seepage flow in a porous medium such as 
soils is based on experiments. A more theoretical approach is to derive this law and predict 
the hydraulic conductivity. This theory is possible if the granular structure is periodic. The 
method of homogenization can then be applied (Sanchez-Palencia, 1974; Benoussan, et. al., 
1978). 

Let a rigid porous medium be saturated by an incompressible Newtonian fluid of constant 
density. Driven by an ambient pressure gradient, the steady flow velocity and pressure p 
in the pores are governed by Navier-Stokes equations 



and 



pu >w j = ~w i ' " ( } 



On the wetted surface T of the solid matrix there is no slip 

Ut = XiET. (6.2.3) 

For slow flows the two terms on the right-hand side of (6.2.2), representing the pressure 
gradient and the viscous force, are equally important. Both the change of pressure and the 
flow velocity vary according to two scales: the microscale t characteristic of the size of pores 
and grains, and the macroscale L imposed by global constraints. Again we assume that 
!/L<Cl. Equating the order of magnitudes of the global pressure gradient P/L to the local 
viscous stress, we get 
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which defines the velocity scale U. Let us normalize the space coordinates by the local scales 
and the unknowns according to the estimates just found 

x i = trf i , p = Pp', Ui = ^«;, (6.2.4) 

where the primed quantities are dimensionless, and p here stands for the change in pressure 
since it only appears in differential form. Equation (6.2.2) becomes, formally, 

where 

o£ 2 P 

Re = e ^_L (6.2.6) 
A* 

is just the Reynolds number. The dimensionless continuity equation remains in the form of 
(6.2.1) and need not be repeated. 

We assume the Reynolds number to be no greater than 0(e), i.e., 

Re < 0(e), or ^ < 0(1). (6.2.7) 
A* 

Note that the pressure gradient term appears formally dominant in (6.2.5), because is 
normalized by the micro-length scale £ in every term. 

From here on we shall return to dimensional variables but retain the ordering symbol e 
to keep track of the bookkeeping 

e V |^ = _|P + eM vV (6.2.8) 

OXj OXi 

The boundary condition on the wetted surface V of the pores is already given by (6.2.3). 

Let us assume that the geometry of the porous matrix is periodic on the microscale, as 
depicted in Figure 6.2.1, although the structure may still change slowly over the macroscale 
L. Each periodic cell Q, is a rectangular box of dimension 0(1). We then expect Ui and p to 
be spatially periodic from cell to cell, while changing slowly over the macroscale. 

We now introduce the multiple-scale coordinates 

Xi, X { = exi (6.2.9) 

and the perturbation expansions 

u i = u? ) +eu\ 1) + (?uf t) + ... (6.2.10) 

p = p(°) + ep ( i ) + e y 2) + ..., (6.2.H) 

where u^\p^ are functions of and 
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Figure 2. A typical O-cell. £2 = Q, U fi 5 . 



Figure 6.2.1: A unit cell on the microscale 
From (6.2.1) we get at the first two orders 0(e°) and O(e) 

(6.2.12) 



dxi 



du {l) duW 
du {2) du { - 1] 



Similarly we get from (6.2.8) 



On the grain/water/interfaces T the velocity vanishes, hence 

uf ] = uf = u 4 (2) = • • • = Xi e r. (6.2.18) 
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In a typical Q cell the flow must be periodic 

uf } uf \ uf\ ■ ■ ■ p (0) p (1) p (1) , • • • are Vt periodic. (6.2.19) 

It may be pointed out that it is the assumption (6.2.7) that renders the perturbations 
equations linear at each order. Had we assumed the Reynolds number to be finite, i.e., 
Re = 0(1), then the convective inertia would be of the order 0(e) and (6.2.16) would be 
replaced by 

Together with (6.2.12) the cell problem is fully nonlinear (Ene and Sanchez-Palencia, 1975). 
From (6.2.15) it is clear that 

p(°)=p(°)(X i ). (6.2.21) 
Because of the linearity of (6.2.12) and (6.2.16) uf and p^ can be formally represented by 

u (0)__ R ,.dP^ p (D = _ A .^ + p(i) f6 2 2 2) 

U, - K tjQX _ p - A Jg ^ +p , yo.l.ll) 

where p^'(Xi) is independent of Xi. It then follows that Kij(xi,Xi) and Aj(xi,Xi) must 
satisfy 

— = (6.2.23) 

OXi 

3 A ■ 

--^- + ^ 2 K ij = -5 ij , (6.2.24) 

where 

Kij = on T (6.2.25) 

Kij, Aj are Q periodic. (6.2.26) 

These four equations define a linear boundary-value problem in an Q cell, which must be 
solved numerically for any prescribed microstructure. 
Defining the average over an cell by 

{f) = W\I fdn > (6 - 2 - 27) 

where flf is the fluid volume inside the $7 cell, we get 

<«?> = -W^ (6-2-28) 



(pW) = -(A J )^ + np^, (6.2.29) 



dXj 
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where n denotes the porosity, which is the ratio of fluid volume in the cell to the total cell 
volume 

Q, 

n=-^-. (6.2.30) 



The Vt average of (6.2.14) gives 



a<«r> i t ^ _ a 



dXi \tt\ J Qf dx 

By virtue of the Gauss theorem and the boundary condition (13.6.30), the volume integral 
vanishes, hence 

0<«i O) > 



dXi 



0, (6.2.31) 



which implies, in turn, that 



-JLflK 
dXi [} ij) 0X 3 ) 



W-^r- =°- ( 6 - 2 - 32 ) 



Equations (6.2.28) and (6.2.31) or (6.2.32) govern the seepage flow in a rigid porous medium 
on the macroscale(Sanchez-Palencia, 1980; Benssousan, et. al., 1978; Keller, 1980). In par- 
ticular (6.2.28) is the classic Darcy's law, where (Kij) is the tensor of hydraulic conductivity. 
If the medium is isotropic and homogeneous on the L scale, we have 

(Kij) = KSij (A,) = 0, (6.2.33) 

where K is a constant. It follows from (6.2.32) that 

ay ) 



dX k 8X k 



0. (6.2.34) 



With proper boundary conditions on the macroscale, p^ can then be found. 

Note that if the Reynolds number is of order unity, the local velocity and pressure p^ 1 
then depend on the global pressure gradient dp(°' /dXi in a nonlinear way through (13.6.32). 
Equations (6.2.22) and Darcy's law (6.2.28) no longer hold. 

In the literature the term homogenization is sometimes used to mean any theoretical 
procedure for deriving macroscale equations by smoothing over the microscale. Homoge- 
nization by multiple scales is just one of the systematic procedures, limited to media with 
a periodic microstructure. In this procedure there is no need for additional closure assump- 
tions. The derivation of effective coefficients is reduced to the solution of certain canonical 
boundary-value problems in a cell. One can also examine the details of the microscale once 
the macroscale problem is solved. 

The cell problem has been solved numerically by finite elements by Sun, Lee and Mei 
(1996, International J. of Heat and Mass Transfer, 39, pp 661-676), using a Wigner-Seitz 
ball in each cell, as shown in Figure 6.2.2. The computed permeability agrees fairly well [ 
with the empirical formula of Carmen and Kozeny based on data for irregular grains. 




Figure 6.2.2: 1.9ft A unit cell having one Wigner-Seitz grain. Csnten 1/Sth of ffie cell with neighboring 
grains in contact overs finite area. Right 1/Sth of the cell with diamond-like grains not in contact, 
i.e.. in suspension. (Adapted from Lee et al. Int. J. Heat 8, Mas Transfer, 33(4), 1956). 



